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We study slightly generalized quantum fidelity susceptibilities where the differential change in the 
fidelity is measured with respect to a different term than the one used for driving the system towards 
a quantum phase transition. As a model system we use the spin-1/2 Ji — J2 antiferromagnetic 
Heisenberg chain. For this model, we study three fidelity susceptibilities, Xp> X-D and xaf, which 
are related to the spin stiffness, the dimer order and antiferromagnetic order, respectively. All 
these ground-state fidelity susceptibilities are sensitive to the phase diagram of the Ji — J2 model. 
We show that they all can accurately identify a quantum critical point in this model occurring at 
J 2 ~ 0.241Ji between a gap less Heisenberg phase for J2 < J 2 and a dimerized phase for J2 > J 2- 
This phase transition, in the Berezinskii-Kosterlitz-Thouless universality class, is controlled by a 
marginal operator and is therefore particularly difficult to observe. 



I. INTRODUCTION 

The study of quantum phase transitions, especially in 
one and two dimensions, is a topic of considerable and 
ongoing interest X Recently the utility of a concept with 
its origin in quantum information, the quantum fidelity 
and the related fidelity susceptibility, was demonstrated 
for the study of quantum phase transitions (QPT) £~— It 
has since then been successfully applied to a great num- 
ber of systems In particular, it has been applied to 
the J\ — J2 model that we consider hereJ^. For a recent 
review of the fidelity approach to quantum phase tran- 
sitions, see Ref. [T3]. Most of these studies consider the 
case where the system undergoes a quantum phase tran- 
sition as a coupling A is varied. The quantum fidelity 
and fidelity susceptibility is then defined with respect to 
the same parameter. Apart from a few studies ^ir— rel- 
atively little attention has been given to the case where 
the quantum fidelity and susceptibility are defined with 
respect to a coupling different than A. Here we consider 
this case in detail for the Ji — J2 model and show that, if 
appropriately defined, these general fidelity susceptibili- 
ties may yield considerable information about quantum 
phase transitions occurring in the system and can be very 
useful in probing for a non-zero order parameter. 

Without loss of generality, the Hamiltonian of any 
many-body system can be written as 



H(X) = H Q + XH X , 



(1) 



where A is a variable which typically parametrizes an in- 
teraction and exhibits a phase transition at some critical 
value A c . In this form H\ is then recognized as a term 
that drives the phase transition,- Using the eigenvectors 
of this Hamiltonian the ground-state (differential) fidelity 
can then be written as: 



f(A) = K*o(A)|*o(A + *A))|. 
A series expansion of the GS fidelity in <5A yields 

(SX) 2 d 2 F 



F(X) = 1 



2 dX 2 



+ 



(2) 



(3) 



where d 2 F = xx is called the fidelity susceptibility. For 
a discussion of sign conventions and a more complete 
derivation see the topical review by Gu, Ref. [l]J If the 
higher-order terms are taken to be negligibly small then 
the fidelity susceptibility is defined as: 



Xx(X) 



2(1 -F (A)) 
(SX) 2 



(4) 



The scaling of xx at a quantum critical point, A c , is often 
of considerable interest and has been studied in detail and 
previous studie d 10 ' 11 ' 14 ! 15 ! 17 have shown that 



Xx~L 2 '\ Xx /N~L 2 i»- 



(5) 



with N = L d the number of sites in the system. An 
easy way to re-derive this result is by envoking finite-size 
scaling. Since 1 — F obviously is dimensionless it follows 
from Eq. Q that the appropriate finite-size scaling form 
for xx is 



XX ~ (6X)- 2 f(L/0. 



(6) 



If we now consider the case where the parameter A drives 
the transition we may at the critical point A c identify 
5X with A — X c . It follows that £ ~ (5X)~ V . As usual, 
we can then replace f(L/^) by an equivalent function 
f(L 1 l v 5X). The requirement that X x remains finite for a 
finite system when SX — > then implies that to leading 



order f(x) 



L 2 / U (SX) 2 , from which Eq. © follows. 



Here we shall consider a slightly more general case 
where the term driving the quantum phase transition is 
not the same as the one with respect to which the fi- 
delity and fidelity susceptibility are defined. That is, one 
considers: 



H(X,5)=H 1 + 5H I , H X =H + XH. 



x- 



(7) 



The fidelity and the related susceptibility is then defined 
as 



F(X,S) = K*o(A,0)|*o(A^))| 
2(1-F(A,<5)) 



s 2 



(8) 
(9) 



2 



The scaling of \8 a t A c for this more general case was 
derived by Venuti et al^ where it was shown that: 



Xs 



L 



2d+2z-2A„ 



Xs/N ~ L 



d+2z~2A v 



(10) 



Here, z is the dynamical exponent, d the dimensionality 
and A„ the scaling dimension of the perturbation Hi. 
In all cases that we consider here z = d = 1. We note 
that Eq. ([TUf assumes [Hi, Hi] ^ 0, if Hi commutes 
with Hi then F = 1 and x<5 — 0- The case where H\ 
and iij coincide is a special case of Eq. (flO|) for which 
Ay- = d + z — 1/f and one recovers Eq. ((5]). 

A particular appealing feature of Eq. ([5]) is that when 
2/v > d, x\/N will diverge at A c and the fidelity sus- 
ceptibility can then be used to locate the A c without any 
need for knowing the order parameter. Secondly, it can 
be shown^ii that the fidelity susceptibility can be ex- 
pressed as the zero-frequency derivative of the dynami- 
cal correlation function of Hi, making it a very sensitive 
probe of the quantum phase transition^ On the other 
hand, if a phase transition is expected one might then 
use the fidelity susceptibility as a very sensitive probe 
of the order parameter through a suitably defined H$ in 
Eq. ([7]). This is the approach we shall take here using 
the Ji — Ji spin chain as our model system. 

The spin- 1/2 Heisenberg Ji — J 2 chain is a very well 
studied model. The Hamiltonian is: 



H 



Si ■ Si+i + J2 Si ■ 



(11) 



where J2 is understood to be the ratio of the next-nearest 
neighbor exchange parameter over the nearest neighbor 
exchange parameter (J2 — Jf/JD- This model is known 
to have a quantum phase transition of the Berezinskii- 
Kosterlitz-Thouless (BKT) universality class occurring 
at Jf between a gapless 'Heisenberg' (Luttinger liquid) 
phase for J2 < Jf and a dimerized phase with a two-fold 
degenerate ground-state for J 2 > Jf . Field theor y 19 ' 20 , 
exact diagonalizatio n 21 i 22 and DMR G 23 i 24 , have yielded 
very accurate estimates of the Luttinger Liquid-Dimer 
phase transition, the most accurate of these being due to 
Eggert 2 ^, which yielded a value of Jf = 0.241167. Pre- 
vious studies by Chen et al*£ of this model using the 
fidelity approach used the same term for the driving and 
perturbing part of the Hamiltonian as in Eq. (jTJ) with the 
correspondence H = J^Si • S i+ i, H\ = J2iSi • S l+2 , 
A = J 2 . 12 . Chen et al. demonstrated that, though 
no useful information about the Luttinger Liquid-Dimer 
phase transition could be obtained directly from the 
ground-state fidelity (and similarly the fidelity suscep- 
tibility), a clear signature of the phase transition was 
present in the fidelity of the first excited stated Some- 
times this is taken as an indication that ground-state 
fidelity susceptibilities are not useful for locating a quan- 
tum phase transition in the BKT universality class. Here 
we show that more general ground-state fidelity suscep- 
tibilities indeed can locate this transition. 

Specifically, we will study three fidelity susceptibilities, 
X p , Xd and xaf, which are coupled to the spin stiffness, 



a staggered interaction term and a staggered field term, 
respectively. In section [H] we present our results for Xp 
while section iHll is focused on xd and section|lV]on xaf- 



II. THE SPIN STIFFNESS FIDELITY 
SUSCEPTIBILITY, Xp 



We begin by considering the Ji — J2 model with J2 = 
but with an anisotropy term A, what is usually called the 
XXZ model: 

ffxxz = £[A£?Sf +1 + \{S+S7 +l + S7S+ + i)]. (12) 



The Heisenberg phase of this model, occurring for A G 
[—1,1], is characterized by a non-zero spin stiffnes s 25 ' 26 
defined as: 



P(L) 



d 2 e{cj)) 



ds 2 



(13) 



6=0 



Here, e(<p) is the ground-state energy per spin of the 
model where a twist of <j) is applied at every bond: 



#xxz(A, < 



E 



[AS?Sf +1 + -(S i S i+l e 



i(/> 



+ SrS+ e" 4 *)]. (14) 



The spin stiffness can be calculated exactly for the XXZ 
model for finite L using the Bethe ansatz;21 and exact ex- 
pressions in the thermodynamic limit are available. 25 i 26 
Interestingly the usual fidelity susceptibility with respect 
to A can also be calculated exactly 28 ' 29 

Since the non-zero spin stiffness defines the gapless 
Heisenberg phase it is therefore of interest to define a 
fidelity susceptibility associated with the stiffness. This 
can be done through the overlap of the ground-state with 
<f) = and a non-zero <fi. With v I / o(A, 4>) the ground-state 
of /Jxxz(A, 4>) we can define the fidelity and fidelity sus- 
ceptibility with respect to the twist in the limit where 
6^0: 



Xp(A) 



K*o(A,O)|* o (A,0))| 
2(1-F(A,£» 



(15) 
(16) 



To calculate Xp the ground-state of the unperturbed 
Hamiltonian was calculated through numerical exact di- 
agonalization. The system was then perturbed by adding 
a twist of e 1 ^ at each bond and recalculating the ground- 
state. From the corresponding fidelity, Xp w& s calculated 
using Eq. (|16|) . Our results for Xp/ L versus A are shown 
in Fig. [TJ For all data </> was taken to be 10~ 3 and pe- 
riodic boundary conditions were assumed. In all cases it 
was verified that the finite value of <j> used had no effect 
on the final results. The numerical diagonalizations were 
done using the Lanczos method as outlined by Lin et al. 30 
Total S z symmetry and parallel programming techniques 




FIG. 1. (Color online.) ^ vs. A: The spin stiffness fidelity 
susceptibility (xp(A)/L) as a function of the z-anistropy A. 
At the A = point the spin-current operator J and kinetic 
energy T commute with the XXZ Hamiltonian and thus such 
a perturbation does not change the ground-state, and the fi- 
delity is one. Thus, Xp m zero a t this point. 



FIG. 2. (Color online.) Xp vs. L (the XXZ model at dif- 
ferent values of A): This graph shows the scaling of Xp with 
system size for different values of the z-anisotropy A. The 
points represent numerical data and the lines represent fits to 
the scaling form predicted for the spin stiffness susceptibility 
Xp/L = A X L + A 2 + A 3 L^ 1 + A 4 L 3 ~ SK . It can be seen that 
there is good agreement. 



were employed to make computations feasible. Numeri- 
cal errors are small and conservatively estimated to be on 
the order of 10" 10 in the computed ground-state energies. 

In order to understand the results in Fig. [T] in more 
detail we expand Eq. (fT4l for small (j>: 



H XXZ (A, 4>) ~ ffxxz(A) + <PJ - 




j= l zY.^ s ^-^ s ti)> (18) 

i 

^E^r+i + sr^+i)- (19) 

i 

Here, J is the spin current and T a kinetic energy term. 
The first thing we note is that, when A = both J and 
T commute with 7Jxxz(A = 0). The ground-state wave- 
function is therefore independent of (for small <j>) and 
X P = 0. This can clearly be seen in Fig. [T] 

In the continuum limit the spin current J can be ex- 
pressed in an effective low energy field theory^ with scal- 
ing dimension Aj = 1. However, we expect sublead- 
ing corrections to arise from the presence of the opera- 
tors (d x $) 2 with scaling dimension 2 and cos(v / 167rK<I>) 
with scaling dimension AK. Here, K is given by K = 
ir/ (2(ir — arcos(A))). For A^O both of these terms will 
be generated by the term T in Eq. (fT7"|). 15 With these 
scaling dimensions and with the use of Eq. (|10[) we then 
find: 

Xp /L = A 1 L + A2 + A 3 L- 1 + A 4 L 3 - SK (20) 

In Fig. [2] a fit to this form is shown for 3 different val- 
ues of A = 0.25,0.5 and 0.75 in all cases do we observe 
an excellent agreement with the expected form with cor- 
rections arising from the last term L 3 ~ SK being almost 



un-noticeable until A approaches 1. We would expect 
the sub-leading corrections L~ x and L 3 ~ 8K to be absent 
if the perturbative term is just <j)J . 

We now turn to a discussion of a definition of xp ln the 
presence of a non-zero J 2 but restricting the discussion 
to the isotropic case A = 1. In this case we define: 

H{4>) = 

E^f+i + \{StSr +l e^ + SrS+ 1 e~ i *)} + 

i 

^ E^+ 2 + \{StS^ + SrS+ie-**)]. (21) 

i 

That is, we simply apply the twist <f> at every bond of the 
Hamiltonian. As before we can expand: 

H(<f>) ~ H(0) + 4>(J X + J 2 ) -t-(Ti +%) + ... ,(22) 

^lE^m-^ii)' ( 23 ) 

i 

z* = \ £^ +5 r +2 _ SrS + a)j (24) 

i 

ri = ^E( 5 ^r+i + ^ + + i)> (25) 

i 

^ = \ E(^ +s r +2 + srs+ +2 ). (26) 

i 

Our results for Xp/L versus J 2 using this definition are 
shown in Fig. [3] for a range of L from 10 to 32. In the 
region of the critical point at J 2 — 0.241167 the size de- 
pendence of Xp/L vanishes yielding near scale invariance. 
How well this works close to J 2 is shown in the inset of 
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FIG. 3. (Color online.) 2j£ vs. J2 and Inset: The generalized 
spin stiffness susceptibility, \ p as a function of the second 
nearest-neighbor exchange coupling J2. The system acquires 
a clearly size invariant form in the vicinity of the critical point 
J2 ~ 0.24 (as well as tending to a global minima). Inset 
shows the the minima for system sizes L=16,20,24,28,32 with 
J| indicated as the vertical dashed line. A clear dependence 
of the J2 value of Xp/L minima on the system size can be 
seen. 



FIG. 4. (Color online.) A: The J2 value of Xp minima as a 
function of system size, as well as a (power law) line of best fit. 
As the system size tends towards infinity the power law best 
fit predicts a minima at J2 = 0.24077 in good agreement with 
previously published results. B: Scaling of Xp at J2 = 0.23 
(the highest, linear curve), J2 = 0.25 (the second highest, 
linear curve) and at the critical point J2 = 0.241167 (flat 
curve). The near constant scaling of 2££ at the critical point 
as well as non-constant scaling on either side of the critical 
point can clearly be seen. 



Fig. [3J This alone can be taken to be a strong indication 
of Xp/ Ls sensitivity to the phase transition. In fact, this 
scale invariance works so well that one can locate the crit- 
ical point to a high precision simply by verifying the scale 
invariance. This is illustrated in Fig. @}3 where Xp/L 1S 
plotted as a function of L for J 2 = 0.23, J2 = Jf and 
J2 = 0.25. From the results in Fig. [4)3 the critical point 
Jf where Xp/L becomes independent of L is immediately 
visible. 

As can be seen in the inset of Fig. [3J Xp/L reaches a 
minimum slightly prior to Jf . The J2 value at which 
this minimum occurs has a clear system size dependence 
which can be fitted to a power-law and extrapolated to 
L = 00 yielding a value of J^ c = 0.24077. Hence, the 
minimum coincides with Jf in the thermodynamic limit. 
This is shown in Fig. Comparison of this value with 
the accepted Jf = 0.241167 reveals impressive agree- 
ment. Another noteworthy feature of the results in Fig. [3] 
is that Xp/L is non-zero at the critical point, Jf. This 
value is very small but we have verified in detail that 
numerically it is non-zero. 

The scale invariance of Xp/L is clearly induced by the 
disappearance^ of the marginal operator cos(v^l67rl?<I>) 
at Jf ■ We expect that in the continuum limit the absence 
of this operator implies that the spin current commutes 
with the Hamiltonian resulting in Xp being effectively 
zero at Jf . The observed non-zero value of Xp/L would 
then arise from short-distance physics. 

Note that, as mentioned previously, we take the spin 
stiffness to be represented by a twist on every bond, both 
first and second nearest neighbor and not merely on the 
boundary as is sometimes done. This choice is not just 



a matter of taste. Imposing a twist only on the bound- 
ary (usually) breaks the translational invariance of the 
ground-state and, through extension, effects the value 
and behavior of the fidelity itself. Another point of note 
is the use of a twist of only <f> between next-nearest neigh- 
bors. Geometric intuition would suggest that a twist 
of 2<j) should be applied between next-nearest neighbor 
bonds. However, for the small system sizes available for 
exact diagonlization it is found that a simple twist of cj) 
on both bonds yields significantly better scaling. 

III. THE DIMER FIDELITY SUSCEPTIBILITY, 

XD 

We now turn to a discussion of a fidelity susceptibility 
associated with the dimer order present in the J± — J% 
model for J2 > Jf . This susceptibility, which we call 
Xd, is coupled to the order parameter of the dimerized 
phase by design. Usually in the fidelity approach to quan- 
tum phase transitions one considers the case where the 
ground-state is unique in the absence of the perturba- 
tion. This is not the case here, leading to a diverging 
Xd/L in the dimerized phase even in the presence of a 
gap. Specifically, we consider a Hamiltonian of the form: 

H = Y}S* ■ Si+i + J 2 Si • s l+2 + sh(-iys t ■ S i+1 ] 

i 

(27) 

Thus, in correspondence with Eq. (J7J we have Hi = 
(—iySi ■ Si+i and we choose the driving coupling to be 




FIG. 5. (Color online.) vs. J2. The generalized 

dimer fidelity susceptibility \d /I/ 3 as a function of the second 
nearest-neighbor exchange parameter J2. A clear intersection 
of all curves can be seen in the vicinity of the proposed crit- 
ical point at J2 ~ 0.2 — 0.25. The inset explicitly shows the 
crossing of L — 12 and L = 24. The dashed vertical lines 
indicate Jf. 

J2. This perturbing Hamiltonian represents a conjugate 
field for the dimer phase. The scaling dimension of Hi is 
known 32 , = |, and from Eq. (fTU|) we therefore find: 

Xd~L 4 - 2A °=L 3 (atJf) (28) 

Due to the presence of the marginal coupling we can- 
not expect this relation to hold for J2 < Jf . However, 
the marginal coupling changes sign at Jf and is there- 
fore absent at Jf where Eq. (|28|) should be exacts For 
J2 < 0.241167 it is known 32 that logarithmic corrections 
arising from the marginal coupling for the small system 
sizes considered here lead to an effective scaling dimen- 
sion Ad > \ • At J2 = Affleck and Bonner— estimated 
Ad = 0.71. Hence, using this results at J 2 = 0, we 
would expect that xd ~ L 2 - 5S which we find is in good 
agreement with our results at J2 =0. 

We now need to consider the case J 2 > 0.241167. At 
J 2 = 1/2 the model is exactly solvable^ 3 , and the two 
dimerized ground-states are exactly degenerate even for 
finite L. For Jf < J2 < 1/2 the system is gapped with a 
unique ground-state but with an exponentially low-lying 
excited state. In the thermodynamic limit the two-fold 
degeneracy of ground-state is recovered, corresponding to 
the degeneracy of the two dimerization patterns. From 
this it follows that xd is formally infinite at J2 = and 
as L —> 00 for Jf < J 2 < 1/2 we expect xd to diverge 
exponentially with L. At Jf we expect xd to exactly 
scale as L 3 and for J2 < J 2 we expect xd ~ L" s with 
a e g < 3. Hence, if xd/L 3 is plotted for different L we 
would expect the curves to cross at Jf. However, the 
crossing might be difficult to observe since it effectively 
arise from logarithmic corrections. 

Our results for xd/L 3 are shown in Fig. [5J where a 
crossing of the curves are visible around J2 ~ 0.2 — 0.25. 
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FIG. 6. (Color online.) A: Scaling of xd vs. L a at the 
critical point J2 = 0.241167. A near perfectly linear scaling is 
observed. B: The J2 value of the intersection of between 
systems of size L and L + 2 plotted as a function of L. The 
curve can be fitted with a power-law line of best fit. The line 
of best fit is found to converge to J2 = 0.241. 

As an illustration, the inset of Fig. [5] shows the crossing 
of L = 12 and L = 24. In order to obtain a more precise 
estimate of Jf the intersection of each curve and the curve 
corresponding to the next largest system were tabulated 
(L and L + 2). These intersection points as a function 
of system size were then plotted Fig. and found to 
obey a power-law of the form a — bL~ a with a ~ 1.8 and 
a = 0.241. This estimate of the critical coupling is in 
good agreement with the value of Jf = 0.241167. 22 

To further verify the scaling of xd at Jf we show in 
Fig. [6j3 xd at Jf as a function of the cubed system size, 
L 3 . The strong linear scaling is in contrast to the scal- 
ing a small distance away from the critical point (not 
shown) where the scaling was found to be distorted by 
logarithmic corrections. 



IV. THE AF FIDELITY SUSCEPTIBILITY, X af 

Finally, we briefly discuss another fidelity susceptibil- 
ity very analogous to xd- We consider a perturbing term 
in the form of a staggered field of the form ^Y^t 
with an associated fidelity susceptibility, xaf- The scal- 
ing dimension of such a staggered field is Aaf = \ and 
as for xd we therefore expect that xaf ~ L 3 at Jf . How- 
ever, in this case it is known^ 2 - that the effective scaling 
dimension for J2 < Jf is smaller than i resulting in 
Xaf ~ L a " !! with a e g > 3 for J 2 < Jf . On the other 
had, in the dimerized phase xaf must clearly go to zero 
exponentially with L. Hence, if xaf is plotted for differ- 
ent L as a function of J2 a crossing of the curves should 
occur. 

Our results are shown in Fig. [7] where xaf/L 3 is plot- 
ted versus J 2 for a number of system sizes. It is clear 
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FIG. 7. (Color online.) xaf/L 3 versus J2. xaf is expected to 
approach zero exponentially with the system size for J2 > J2, 
to scale as L 3 at Jf and to scale as L a <= a with ot e g > 3 for 
J2 < J%- A crossing close to the critical point J 2 (dashed 
vertical line) is then visible. 



ity beyond a simple perturbation of the same term that 
drives the quantum phase transition. By using the spin- 
1/2 Heisenberg spin chain as an example we first cre- 
ated a susceptibility which was directly coupled to the 
spin stiffness but of increased sensitivity. This fidelity 
susceptibility, which we labelled Xp can be used to suc- 
cessfully estimate the transition point at J2 ~ 0.241. 
Next we constructed another fidelity susceptibility, xd, 
this time coupled to the order parameter susceptibility 
of the dimer phase. Again, we were able to estimate the 
critical point at a value of 0.241. Finally, we discussed 
an anti ferromagnetic fidelity susceptibility that rapidly 
approaches zero in the dimerized phase but diverges in 
the Heisenberg phase. Although susceptibilities linked to 
these quantities appeared the most useful for the J\ — J2 
model we considered here, it is possible to define many 
other fidelity susceptibilities that could provide valuable 
insights into the ordering occurring in the system being 
studied. 



from these results that xaf indeed goes to zero rapidly 
in the dimerized phase as one would expect. Close to 
J| the scaling is close to 1? where as for J2 < Jf it is 
faster than L 3 . Hence, as can be seen in Fig. [JJ a crossing 
occurs close to J$. 

V. CONCLUSION AND SUMMARY 

In this paper we have demonstrated the potential ben- 
efits of extending the concept of a fidelity susceptibil- 
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